The statistics of Wigner delay time in Anderson disordered systems 
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We numerically investigate the statistical properties of Wigner delay time in Anderson disordered ID, 2D and 
quantum dot (QD) systems. The distribution of proper delay time for each conducting channel is found to be 
universal in 2D and QD systems for all Dyson's symmetry classes and shows a piece-wise power law behavior 
in the strong localized regime. Two power law behaviors were identified with asymptotical scaling r~^'^ and 
T~'^ , respectively that are independent of the number of conducting channels and Dyson's symmetry class. Two 
power-law regimes are separated by the relevant time scale tq ~ ft/ A where A is the average level spacing. It 
is found that the existence of necklace states is responsible for the second power-law behavior t^^, which has 
an extremely small distribution probability. 
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I. INTRODUCTION 

After the pioneering worksri^ the problem of quantum me- 
chanical scattering has attracted intensive research interest in 
many fields. The dynamical aspect of the scattering process 
can be characterized by the energy derivatives of the scattering 
matrix, known as the Wigner delay time r. Semi-classically, 
Wigner delay time measures the time spent by the center of 
a wave packet in the scattering region, which is simply re- 
lated to its group velocity. Since r is critically dependent on 
the transport process, it is not self-averagingi^. Hence a com- 
plete distribution is required to fully comprehend this quantity 
in classical chaotic or quantum disordered systems. 

The distribution of delay time in ID system has been thor- 
oughly studied both theoreticallyi^"— and experimentally.iJ^^ 
Texier and Comet'- showed by various methods the universal- 
ity of the distribution of r in a ID semi-infinite system, which 
has a r^^ power-law tail at large t in the localized regime. 
This power law tail has also been confirmed by others ji^iiii^ 
It has been established that r^^ behavior is valid for different 
types of disorder potential with 5, gaussian, box, and expo- 
nential decaying distributions. ""'- 

For 2D systems with multichannels, most of the theoretical 
works are based on the random matrix theory (RMT)ifi2r2^ A 
simple expression of Pjv {t) for different Dyson's symmetry 
classes has been derived-^'?^ at the ideal coupling condition, 
which has a distribution with I3{N + 1) degrees of free- 
dom where /3 is the symmetry index and N the number of 
conducting channels. On the other hand, a P{t) ^ r^'^-^ 
asymptotic behavior in the diffusive regime was numerically 
obtained from the kicked rotor model (KR model)i21— Start- 
ing from RMT, it was shown that an open chaotic system 
weakly coupled with many open channels also exhibits the 
T^^-^^ power-law behavior in symmetry classes /3=1, 2, and 
4^25j26 5q f^j. ^jj jjjg studies are focused on diffusive regime, 
less attention has been paid on the localized regime. It is the 
purpose of this paper to fill this gap. In this paper, we wish to 
explore the possibility of universal behavior of the distribution 
of Wigner delay time r in the localized regime. 

In this paper, we carry out an extensive numerical inves- 
tigation of statistical properties of Wigner delay time r in 
disordered systems with Anderson-type of impurities for ID, 
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FIG. 1: sketch of the three geometries of interest, (a) ID semi-infinite 
tight-binding chain; (b) 2D system; (c) QD system. 



2D and quantum dot (QD) systems. For ID systems, our re- 
sults confirm that the distribution of Wigner delay time fol- 
lows a power-law behavior r^^ in the localized regime. For 
2D and QD systems, our results show that the distribution of 
proper delay time of each conducting channel obeys a uni- 
versal piece-wise power-law in strong localized regime that 
is independent of number of channels and Dyson's symme- 
try index j3. For a scattering system, the characteristic time 
scale To is related to the group velocity Vg of the electron, i.e., 
Tq ^ L /vg where L is the characteristic length of the system. 
Our results show that when r < tq, the distribution follows 
a power-law r^^ '"' while for r > tq a new power-law of 
is obtained for the proper delay time. Our result indicates that 
the power-law of t^^ can only be observed in the localized 
regime. The physical origin of the new power-law behavior is 
the existence of the so called Azbel resonant state- or neck- 
lace state^^i^. This necklace state has a very long lifetime and 
is the multi-resonant state inside the scattering region in the 
localized regime. When the incoming electron has N con- 
ducting channels, the distribution for total delay time is given 
by P{J2i '''i)- Although P{Ti) obeys the piece-wise power- 
law behavior for each proper delay time t^, the distribution 
for total delay time behaves differently. Our results show that 
as N increases the power-law region for t^^-^ becomes nar- 
row and eventually diminishes in the large N limit while the 
power-law region for r^^ remains but shifts towards small t. 
Therefore, in the large N limit, only one power-law scaling of 
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r ^ survives in the localized regime. This conclusion is valid 
for 2D and QD systems and for three different ensembles with 

/3=1,2,4. 

This paper is organized as follows. In section II, a theo- 
retical formalism and a numerical implementation scheme for 
calculating Wigner delay time are given. In section III, exten- 
sive numerical results and analysis for distribution of Wigner 
delay time are presented for ID, 2D and QD systems with dif- 
ferent symmetries {/3 = 1,2,4). Finally a brief summary is 
given in section IV. 



gauge with B the magnetic field, cr is the Pauli matrix and 
tgo is the strength of spin-orbit coupling, fj, = gfis/'^ is the 
magnetic moment, with (7 = 4 the Lande g-factor and ^ b the 
Bohr Magneton. The tight-binding Hamiltonian with nearest- 
neighbor hopping has an expression^ 
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n. THEORETICAL FORMALISM 

Following Wigner- and Smith,"* the Wigner-Smith delay 
time matrix is defined in terms of the scattering matrix S as 



Q{E) ^ ~ihS\E) 



dSjE) 
dE 



The delay time t is simply the summation of the diagonal 
elements of matrix Q(E) 



r = Tr[Q]^-^hTr[S^—] 



(1) 



From now on we will omit the energy dependence of the rele- 
vant quantities for simplicity. Supposing there are N conduct- 
ing channels, the eigenvalues of delay time matrix Q are called 
the proper delay times, ri , T2 , . . . , tat, which can be viewed as 
the contribution to the total delay time r from the correspond- 
ing conducting channel. 

Three setups under investigation are schematically illus- 
trated in FiglH which are respectively (a) a ID semi-infinite 
chain with length L = 1000a lattice points; (b) a 2D sys- 
tem with 80 by 80 lattice sites connected to a single lead with 
width Wo = 80a; (c) a quantum dot system with the same 
number of lattice points as that of the 2D system but with 
a lead of narrower width Wq = 20a. Here a is the lattice 
spacing between two adjacent sites, being the length scale in 
the calculation. Most of numerical calculations were done us- 
ing these parameters. In the present work we use the con- 
ventional nearest-neighbor tight-binding approximation and 
Green's function formalism to numerically study these sys- 
tems. All these geometries are connected by only one semi- 
infinite lead to the electron reservoir, which ensures the uni- 
tary of the scattering matrix since all electrons incident will 
be reflected back into the reservoir in these one-lead systems. 

Assuming our 2D system is in the x-y plane. In the pres- 
ence of perpendicular magnetic field and Rashba spin-orbit 
interaction, the generalized Hamiltonian is given by 



H 



:^[p+-r + V + fiB-a + ^[a X (p + ^)].(2) 
2m* c he 



where p is the momentum and m* the effective mass of elec- 
tron. Here V is the confining potential that is set to zero inside 
the device and infinity at the boundary of the device except at 
the interface of the lead. The vector potential A due to the 
magnetic field is expressed as A=(-By, 0, 0) under Landau 



where c\^^^{cnma) is the creation (annihilation) operator for 
an electron on lattice site (n, m). represents the on-site 
energy, with magnitude 2t for ID chain and 4t for 2D square 
lattice. Here t = h'^/2m*a? is the nearest-neighbor hopping 
energy, which is the energy scale in this work. Disorder en- 
ergy e„„i, which is Anderson-type with a uniform distribution 
in the interval [-W/2,W/2], is added to the on-site energy 
with W the disorder strength. 

Based on the tight-binding Hamiltonian, the retarded 
Green's function in real space is defined as 

C'iE) = {E-H-Y.'')-^ 

where S'' is the self-energy of the lead, which can be calcu- 
lated by a transfer matrix method^. E is the electron Fermi 
energy and chosen to be at the center of the corresponding 
subband. By the Fisher-Lee relation- which connects the scat- 
tering matrix and the Green's function, the delay time r is 
rewritten in terms of G"" as 



oE 



(3) 



where G" is the advanced Green's function, G° = (G'')^, and 
r is the line width function describing coupling of the lead 
to the scattering region which is given by F = — E"]. 
In the presence of large disorders, the system can either be in 
diffusive or localized regimes depending on the dimensionless 
localization length £,/L. In the localized regime, the localiza- 
tion length can be defined as^ < G >= Gcxp(— 2L/^) with 
G the conductance and C is a constant to be determined. To 
eliminate C, the localization length ^ at a particular disorder 
strength can be obtained from 
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In <Gi> -In < G2 > 



From Eq.([3j, it is clear that the calculation of Wigner delay 
time is equivalent to that of density of states. In our numerical 
calculation, the real space Green's function can be obtained 
by matrix inversion which is very time-consuming. For delay 
time, one only needs the first Niead columns of Green's func- 
tion where Niead is the dimension of the line width function 
r. The transfer matrix method is suitable for this purpose and 
is fast. To study the statistics of r, we need to generate an 
ensemble of different realizations of the disordered systems. 
As we will show in the next section, the distribution of t has a 
new power-law for large delay time r > tq. These large delay 
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times correspond to rare events with extremely small proba- 
bility, which means that to study the new power-law regime 
of P(t) a large configuration ensemble is necessary. There- 
fore intensive computation is required to accumulate enough 
data for statistical analysis. To speed up the calculation, we 
can rewrite the line width function F as^ 

N 

r = i[j:- -j:-]=J2\w,){w,\ 

i 

where N is the number of conducting channels in the lead and 
\Wi > is the renormalized eigenfunction of F. Then the delay 
time can be expressed as 

N 

T = hTriGTG''] ^ h^{G''\W^)){G''\W^)y (4) 

i 

This representation can speed up the numerical calculation 
and also enables to calculate the proper delay time from the 
i — th conducting channel. Our results show that using the LU 
decomposition with multi-frontal algorithm to solve the linear 
equation {E — H — = \ Wi > is faster than the transfer 

matrix method especially for large system size like 100 x 100. 
So far we have discussed the algorithm to study the statistics 
of r. The numerical results and relevant discussions is the 
content of the next section. 



ni. NUMERICAL RESULTS AND DISCUSSION 

In this section, we will present our numerical results. To test 
our code, we first show results of ID systems with orthogonal 
symmetry, in which case the distribution of r has been stud- 
ied thoroughly. Then in the following two subsections, we 
shall discuss in detail the cases of 2D and QD systems with 
different symmetries and number of conducting channels. 



A. ID tight-binding cliain 

For the ID system, Texier and Comeii^ derived an analytic 
expression of the distribution of Wigner delay time in the lo- 
calized regime for high energies or weak disorder strengths. 
It was found that P{t) has an algebraic tail in the localized 

29 

regimoai 

P(r) = -^e-^/^^ (5) 

where ^ is the localization length and v is the group velocity. 
They also found numerically that in the ballistic regime P(r) 
obeys a Gaussian distribution. The disorder used in their nu- 
merical simulation is the S potential with random position in 
the system. 

We have calculated the delay time distribution for a ID 
tight-binding chain in both the ballistic regime and the local- 
ized regime at a relative high energy _E = 1.0 with disorder 
strength W — 0.01 and W — 1.1, respectively. The result is 
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FIG. 2: The distribution of Wigner delay time of a ID semi-infinite 
long chain at disorder strength W — 1.1. Inset is the normalized 
distribution histogram of r at = 0.01. Blue curve shows the 
Gaussian fitting of the histogram. 4,000,000 configurations have 
been generated at each disorder strength W. 

shown in Fig[2l Clearly the delay time distribution has a Gaus- 
sian shape at weak disorder W = 0.01 when L ^ ^ (see inset 
of Fig|2]i. As W increases, P{t) transforms gradually from a 
symmetric Gaussian distribution to a one-sided distribution^. 
To make the behavior more transparent, one may change the 
variable from t to its natural logarithm. Take Eq.© as an ex- 
ample. When we use Inr as the variable, P(lnT) ~ tP{t), 
and taking logarithm at both sides, one arrives at 

£ £ 1 

InP(lnT) =ln- -Inr- -- (6) 

V V T 

The first term in the above expression is a constant for a 
specific disorder strength and the last term tends to zero 
when T is very large. As a result, a linear tail arises in the 
In P(ln(r)) — In r curve. In Fig|2]we see that the linear tail of 
lnP(ln(r)) at disorder strength W ~ 1.1, where the system 
is already localized (L ^ £), can be well fitted by a straight 
line with a slope —1, which implies that P{t) oc 1/t^ in the 
large t region. We notice that there is a departure from the 
straight line at the end of the distribution curve. This depar- 
ture can be improved by using more configurations, since the 
larger the t, the smaller its occurrence probability. We will 
explain the origin of these large r in the next subsection. 

To summarize briefly, our numerical results for Anderson- 
disordered ID system are consistent with the conclusion of 
Texier and Comet which confirms the universality of the 
1/t^ power-law tail. 

B. 2D square lattice 

For 2D systems, most theoretical works of delay time dis- 
tribution were within the random matrix theory. As we know, 
the random matrix theory works well in the diffusive regime. 
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FIG. 3: Panel (a) is the distribution histogram of r at disorder strength W = 0.0001 and the blue curve is a standard Gaussian fitting. Panel (b) 
is P{r) atW = 15 and the corresponding histogram with variable changed to In t is displayed in panel (c). Panel (d) contains the In P(ln r) 
versus In r curves for W = 5, 15, 25, and 40. Panel (e) and (f) are respectively the asymptotic behavior of slopes for different linear sections 
in panel (d)'s InP(lnT) curves. The system is under such conditions: N=l and 13 = 1. 



However, it becomes difficult to describe properties in the lo- 
calized regime. In the diffusive regime, RMT predicts that 
P{t) shows2Si2^ a universal power-law behavior r^^-^, which 
was confirmed by a KR model2i2?. Our numerical results for 
proper delay time obtained from 2D Anderson systems also 
show the r^^^ power-law behavior in the intermediate range 
of T and reveal new properties of the distribution of Wigner 
delay time in the localized regime. 

First we start with the simplest case of single conduct- 
ing channel iV = 1 with preserved time-reversal symmetry 
{/3 ~ 1) and the results are shown in Figl3] Fermi energy 
of electron is chosen to be at the center of the first subband, 
E = 0.004. From panel (a) of Fig|3]it is obvious that in the 
ballistic regime the distribution of Wigner delay time t has 
a natural Gaussian shape at a weak disorder W = 0.0001, 
which is the same as that of ID systems. As the disorder 
increases, P{t) is no longer symmetrically distributed but 
spreads over a wide range with a one-sided peak located at 
small delay time region, as shown in panel (b) with W = 15. 
From the calculation of localization length ^, we know that 
the system is localized at such a disorder strength. Therefore 
to get an overview of the distribution including the long delay 
time tail, we change the variable from r to In t and the his- 
togram of P(lnT) is depicted in panel (c). Based on Eq.® 
and the argument therein, we plot the logarithm of the dis- 
tribution of In T in panel (d) for different disorder strengths 
ranging from W = 5toW — 40. We see that as the disorder 



is increased in the scattering region, the curve of InP(lnT) 
vs In r gradually evolves into a piece-wise power-law pattern 
with two different power-laws and an abrupt change from one 
to the other at a particular value of tq. This picture becomes 
clear when W increases to 40. We understand this behavior as 
follows. In the scattering system, the characteristic time scale 
is set by To = L/vg with L 3Lq where Lq is the dimension 
of the scattering region and Vg is the group velocity. In the 
strong localized regime, there are two kinds of scattering event 
naturally separated by tq, one coiTesponds to the usual direct 
reflection with delay time t < tq and the other corresponds 
to multi-resonant reflection with long delay time t > tq. For 
Figl3]with Lq = 80, tq = 3Lo/vg gives Iiitq ~' 7.5 which 
is close to the transition point from the power-law t^^-^ to 
the power-law r^^ in Fig. 3. For strong disordered systems, 
the large delay time scattering events with r > tq are really 
rare events with extremely small probability. In addition, this 
probability decreases as W increases. For instance, the prob- 
ability P{T)\r>ro is 13% for W = 15 while for W = 40 it 
drops to 0.32%. To study such rare events, more and more 
configurations of disordered samples are required for statistic 
analysis to get an accurate result. For = 40 we have used 
an ensemble of over 40,000,000 different realizations. 

Another point worth noticing is that the exponents of two 
power-laws decrease slowly as disorder strength is increased. 
These two exponents versus disorder strength W are plotted in 
panels (e) and (f), which converge towards -0.5 and -1.0, re- 
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FIG. 4; Statistics of Wigner delay time at disorder strength W = 100 
for N — 1. Panel (a) is the P(lnr) histogram and panel (b) is the 
corresponding In P(ln r ) — In r curve. Red and purple straight lines 
individually have slope -0.5 and -1.0. 



spectively. We note that the curve in panel (f) has larger error 
bars compared with that in panel (e). This is understandable 
since large t is more difficult to sample. 

Since the computation is getting extremely time-consuming 
it is very difficult if not impossible to obtain the plot shown in 
panels (e) and (f) for large disorder strengths. Hence we de- 
cide to calculate P(r) for a particular disorder strength which 
is large enough to determine the exponents for both power-law 
behaviors. After some trial and error, we found the disorder 
strength W = 100 is appropriate for this purpose. The expo- 
nent of power-law in the intermediate region is already con- 
verged at such a W and does not change upon further increas- 
ing of disorder strength to = 200. To get a clear view of the 
tail with small probability, a large ensemble of 420,000,000 
configurations is accumulated and the analyzed data is shown 
in FiglD Clearly the algebraic tail where t > tq is invisible in 
the histogram of panel (a), since it accounts for only 0. 1 18 % 
in the ensemble. From panel (b) one clearly sees that the first 
power-law of In P(ln t) has a slope of -0.5, corresponding to 
a power law distribution P(r) r~^-^. After a transition 
region around t = tq shown by the blue circle in Fig HI the 
second power-law distribution is found to have a slope of -1, 
which means that it is a power-law tail P{t) ~ t~^. The 
fluctuation at the end of the tail arises because there is not 
enough configurations. In such a strong localized regime with 
W = 100, delay time with magnitude lur > 15 often has a 
few configurations out of the total 420,000,000 ensemble. 

The power-law distribution r^^^ has already been pre- 
dicted by RMT^ and was also confirmed numerically us- 
ing KR model22 in the diffusive regime. Our results show that 
this power law t~^-^ for each proper delay time exists for 
any where A^ is the number of conducting channels (also 



see numerical results presented below). In addition, this be- 
havior persists in the strong localized regime. Importantly, 
the T^^ algebraic tail is a novel feature that has not been re- 
ported before in 2D systems. The physical origin of the novel 
power-law behavior is due to the Azbel resonance states or 
necklace states^—, which are extended states in localized 
systems. Since in the localized regime most of the electron 
states are localized, these necklace states or Azbel resonances 
survive through multiple-resonances. These states are really 
rare events in the disordered samples so that they have an ex- 
tremely small probability to occur. For example, in the above 
calculation with N=l and W = 100, the probability for those 
states with r > exp(8) is 0.03%. Despite of their rare na- 
ture, these states can play a significant role in the distribution 
function of delay time and dominate in the density of states 
in disordered systems. It is because a necklace state experi- 
ences multiple resonant scattering giving rise to a large den- 
sity of states. For instance, emittance calculations in disor- 
dered ID, 2D, and QD systems show that due to the necklace 
states, the average emittance remains negative in the localized 
regime;;^ In fact, the existence of necklace states has already 
been observed through optical experimental in a quasi- 1 D sys- 
tem. The fact that the signature of necklace states is observed 
again in disordered 2D systems with single lead indicates that 
the non-localized necklace states are generic in strong disor- 
dered systems, although difficult to see them. 
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FIG. 5: Distributions of proper delay time n and T2 and the total r 
at disorder strength W = 100 in the case of 2 conducting channels. 
Inset panel (a) highlights the tirst linear section of distribution with a 
-0.5 slope (blue-dash straight line), while panel (b) for the algebraic 
tail with a -1 slope (purple-dash-dot-dash straight line). Ensemble 
size is 150,000,000. 

Next we consider the case of two conducting channels 
(A^ = 2) in a system with (3^1. Here electron Fermi energy 
is fixed at the center of the second subband, where E = 0.01. 
Since the group velocity of electron incident from different 
subbands are different, proper delay times for different con- 
ducting channels are not statistically equivalent. Hence the 
properties of the proper delay time t;, which describes the 
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scattering of electron of the i channel, as well as the total 
Wigner delay time r = X^i^'''* ^re studied separately. To 
save computational time, we calculated the distribution of rel- 
evant times only at a large disorder strength W ~ 100. As we 
have seen above, at this W P{t) of N = 1 converges to the 
piece-wise power-law (t^^-^ and r^^) behavior. Statistically 
analyzed result upon an ensemble of 150,000,000 configura- 
tions is shown in Fig|5] To separate the curves we zoom in 
the picture around two particular areas of the two power-law 
regions, which are highhghted in the insets. 
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FIG. 6: Panel (a) shows the distribution of proper delay time ri, rs, 
and Tio as well as the total delay time t dXW = 200 with 10 con- 
ducting channels in the lead. Panel (b) contains the normalized dis- 
tribution histogram and its logarithm correspondence of total delay 
time T at W = 50 for AT = 35. 

From Fig|5]it is clear that proper delay time ti and T2 fol- 
low the same distribution with different constants (P(ti/2) ~ 
Ci/2T°' with a = -1.5 or -2). It shows that the proper delay 
time of the first and second subbands are statistically inde- 
pendent. Since the total delay time r is the summation of ti 
and T2, it is not surprising that the distribution of total delay 
time r keeps the essential feature of individual proper delay 
time. Clearly all three curves show piece-wise power-law be- 
haviors, which is similar to the case of a single channel = 1 
(seeFig|4]i. In view of the distribution for systems with = 1 
and = 2, we expect that the delay time distribution in mul- 
tichannel case {N > 2) also show piece-wise power-law be- 
havior. However the transition region in the distribution curve 
of T is broadened due to the overlap of P{ti) and P{t2). As 
the channel number increases, this may lead to a significant 
change of the distribution form as will be discussed in detail 
below. 



Numerical results for cases of more conducting channels 
iV = 10 and 35 within orthogonal ensemble (/3 = 1) are 
shown in Fig|6] Following observations are in order (1). For 
a strong disorder strength W=200, statistical analysis on the 
distribution of proper delay time at = 10 shows that the 
individual Ti exhibits a piece-wise power-law behavior with 
two scalings and t^^, similar to that of the proper de- 

lay time in A^ = 2 case. In addition, the total delay time 
distribution P{t) ai N = 10 given by the superposition of 
all the P{Ti) (i=l,2,.., 10) results in the same behavior (see 
panel (a)). We also note that due to the group velocity differ- 
ence for different subbands, the transition point tq = 3io / Vg 
of proper delay time distribution is not at the same position 
in liiP(liiT)-lnT curve. The transition value of proper de- 
lay time is the smallest for the first subband while the largest 
for the highest subband. The overlap of proper delay time 
distribution causes the broadening of transition region of the 
distribution P{t). (2). Comparing the situations of Af = 2 
and A^ = 10 (Fig|5]and panel (a) of Fig|6|l, the broadening 
is more obvious at large N. When the number of conducting 
channel increases to A^ = 35, the broadening is so signifi- 
cant that it destroys the power-law behavior of r^^^ in the 
total delay time distribution although t^^ '^ power-law exists 
for each proper delay time distribution. Therefore our results 
suggest that in the large N limit, only the power-law of r^^ 
survives. 

Now we turn to the systems with symmetry class f3 = 2 
where the time reversal symmetry is broken and /? = 4 in the 
absence of spin-rotation symmetry. In our numerical calcula- 
tion, the unitary symmetry class {(3 ~ 2) is realized by apply- 
ing an external magnetic field while for symplectic symmetry 
(/3 = 4) we consider the spin-orbit interaction. In both uni- 
tary and symplectic cases we mainly numerically investigate 
the statistics of Wigner delay time for single conducting chan- 
nel (A^ = 1) since the computation is getting extremely time- 
consuming for 2D with /3 ~ 2 and 4. From the discussion 
in orthogonal symmetry case, some intuition can be obtained 
for the distributions of delay time at multi-channel case. From 
Fig|2]we see that there are clearly two power-law regions with 
exponents approaching —1.5 and —2 for total delay time dis- 
tribution when N=l (panels (a) and (b)). We have tested two 
different magnetic field strengths and they give similar results. 
In addition, the transition phenomena from one power-law to 
another can also be seen from inset of Fig|7ja) and (b) that 
are different from the case with /3 = 1. We have checked 
that the proper delay time distribution follows a piece-wise 
power-law behavior. From our experience in orthogonal en- 
semble we expect that the r^^^ intermediate region will be 
destroyed at large N. Indeed, our results confirm this expecta- 
tion (see panel (c) where A^ = 20 is investigated). Except for 
this difference shown by green circles in Fig|7] the delay time 
distribution is similar to that of the system with reserved time- 
reversal symmetry (/3 = l)(panel (a) of FiglU. Finally, in the 
case of symplectic symmetry, the distribution is the similar to 
that of /3=1. We have also varied the strength of spin-orbit 
coupling tso and found that it does not affect the distribution 
behavior 

To briefly summarize our results for 2D systems, the 
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FIG. 7: Panel (a) and (b) are respectively the distributions of delay time at two different sets of parameter. In panel (a) it is VK = 100 and 
B — 0.01, while in (b) W — 400 and B — 0.02. Both systems have single conducting channel(A'' = 1) in unitary symmetry case (/3 = 2). 
Ensemble size is 63,000,000. Panel (c) shows the statistics of r at multi-channel case with N — 20 al l3 = 2. Case of symplectic symmetry 
case (/? — 4) is drawn in panel (d) with W = lOOand A'^ = 1, where an ensemble of size 120,000,000 configurations is used. 



Wigner delay time distribution in the 2D lattice system shows 
a piece-wise power-law feature in strong localized regime, 
with power-law converging to t~^-^ in the intermediate re- 
gion and the algebraic tail obeying a different power law t'^ 
for a single conducting channel. The two power-law regions 
are separated by a clear transition region determined by tq in 
the distribution curve. For multi-conducting channels, the dis- 
tribution of each proper delay time is found to have the same 
power-law distribution as that of the delay time r for = 1 . 
For a few conducting channels, the distribution of total Wigner 
delay time behaves like that of a single channel. When N is 
large, however, the superposition of each proper delay time 
distribution gives rise to a broad transition region that gradu- 
ally destroys the r^^-^ power-law region. In the large N limit, 
only one power-law region exists which corresponds to r^^. 
These features seem to be independent of symmetry class of 
2D systems, i.e., regardless of orthogonal, unitary or symplec- 
tic symmetries. 



C. Quantum dot system 



The quantum dot (QD) system consists of a square scatter- 
ing region and single lead attached at the right side with 1/4 
width of the scattering region, as shown in Fig[T] Numerical 
results depicted in Fig. 8 show that the conclusions obtained 
in 2D system are also applicable to QD case. In FigUJa) and 
(b), we see that the piece-wise power-law behavior is observed 
again in all Dyson's symmetry classes for the case of = 1. 
For /3 = 1 and N = 18, we see that the power-law r~^-^ 
no longer exists that is similar to the situation in 2D systems. 
A slight difference compared to the 2D case is that there is 
an additional peak in the distribution of delay time P(t) vs 
T that indicates the onset of the power-law region of T^^(see 
panel (c) of Fig [8]). To check the size effect of the distribution, 
we have examined a larger quantum dot system with size of 
L = 100 and lead width Wq = 10. The numerical results for 
f3 — 1 and = 1 are shown in Fig. 8(d) and similar conclu- 
sion can be drawn. 
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FIG. 8: Quantum dot case. Panel (a) and (b) are the distribution of delay time for different symmetry classes (/3 = 1, 2, and 4) with single 
conducting channeKA^ — 1). Panel (c) corresponds to A'^ = 18 in orthogonal symmetry /3 = 1. Calculation parameters are shown in the 
graphs. The result of a large QD at /3=1 is presented in panel (d). 



IV. CONCLUSION 

In conclusion, we numerically investigate the statistical 
properties of Wigner delay time in Anderson disordered ID, 
2D and quantum dot systems for different symmetry classes 
/3 = 1,2,4. The proper delay time distribution is found to be 
universal in strong localized regime for 2D and QD systems 
and shows a piece-wise power-law behavior. In addition to the 
known power-law scaling region r^^-^, a novel power-law re- 
gion of is identified in the localized regime, which is inde- 
pendent of the number of conducting channel N and Dyson's 
symmetry classes /3. Our results indicate that the existence of 
necklace states is responsible for the algebraic tail, which 
are rare events and have an extremely small distribution proba- 



bility. The total delay time distribution can behave differently. 
For a few conducting channel, there is a crossover region from 
one power-law region to another. As the number of conduct- 
ing channels N increases, this crossover region broadens and 
the power-law region of r"^ '"' becomes narrow. In the large 
N limit, the power-law of r~^-^ is destroyed due to the broad- 
ening of the crossover region and only the t^^ algebraic tail 
survives. 

V. ACKNOWLEDGMENTS 

This work is supported by RGC grant (HKU 705409P) from 
the HKSAR and LuXin Energy Group. We thank the HPC 
POWER of the computer center, HKU for computation re- 
sources allocation. 



* Electronic address: jianwang@hkusub.hku.hk 

' L. Eisenbud, Ph. D dissertation, Princeton University, 1948. 

2 E. P Wigner, Phys. Rev. 98, 145 (1955). 

^ R T. Smith, Phys. Rev. 118, 349 (1960). 

"* A. Z. Genack, P. Sebbah, M. Stoytchev, and B. A. van Tiggelen, 



Phys. Rev. Lett. 82, 715 (1999). 
^ B. A. van Tiggelen, P. Sebbah, M. Stoytchev, and A. Z. Genack, 

Phys. Rev. E. 59, 7166 (1999). 
" Z. Qiao, Y. Xing, and J. Wang, Phys. Rev. B. 81, 0851 14 (2010). 
^ M. P. Lopez Sancho, J. M. Lopez Sancho and J. Rubio, J. Phys. 



9 



F: Met. Phys. 14, 1205 (1984); M. P. Lopez Sancho, J. M. Lopez 
Sancho and J. Rubio, J. Phys. F: Met. Phys. 15, 851 (1985). 
D. S. Fisher, and P A. Lee, Phys. Rev. B 23, 6851 (1981). 
' A. MacKinnon and B. Kramer, Phys. Rev. L 47, 1546 (1981); A. 
MacKinnon and B. Kramer, Z. Phys. B: Condens. Matter 53, 1 
(1983). 

'° J. Wang, and H. Guo, Phys. Rev. B 79, 0451 19 (2009). 
" A. Comet and C. Texier, J. Phys. A 30, 8017 (1997). 

S. K. Joshi, A. K. Gupta and A. M. Jayannavar, Phys. Rev. B 58, 

1092 (1998). 

" C. Texier and A. Comet, Phys. Rev. Lett. 82, 4220 (1999). 

F. Steinbach, A. Ossipov, T. Kottos, and T. Geisel, Phys. Rev. Lett. 
85,4426(1999). 

A. Ossipov, T. Kottos and T. Geisel, Phys. Rev. B 61, 
11411(2000). 

S. A. Ramakrishna and N. Kumar, Eur. Phys. J. B 23, 515(2001). 

" W. O. Amrein, and Ph. Jacquet, Phys. Rev. A 75, 022106(2007). 
J. D. Bodyfelt, J. A. Mendez-Bermiidez, A. Chabanov and T. Kot- 
tos, Phys. Rev. B 77, 045103(2008). 

" A. Z. Genack, P. Sebbah, M. Stoytchev, and B. A. van Tiggelen, 
Phys. Rev. Lett. 82, 715 (1999). 

^° A. A. Chabanov and A. Z. Genack, Rev. Lett. 87, 233903 (2001). 

^' K. A. Muttalib and P Wolfle, Phys. Rev. Lett. 83, 3013 (1999). 
Y. V. Fyodorov and H.-J. Sommers, Phys. Rev. Lett 76, 4709 
(1996). 



P Seba, K. Zyczkowski and J. Zakrzewski, Phys. Rev. E 54, 2438 
(1996). 

P W. Brouwer, K. M. Frahm, and C. W. J. Beenakker, Phys. Rev. 
Lett 78, 4737 (1997). 

Y. V. Fyodorov, D. V. Savin and H.-J. Sommers, Phys. Rev. E 55, 
R4857 (1997). 

H.-J. Sommers, D. V. Savin, and V. V. Sokolov, Phys. Rev. Lett 

87, 094101 (2001). 
" A. Ossipov, T. Kottos and T. Geisel, Europhys. Lett 62, 719(2003). 

T. Kottos, A. Ossipov and T. Geisel, Phys. Rev. B 68, 

066215(2003). 
2' T. Kottos, J. Phys. A: Math. Gen. 38, 10761(2005). 

A. D. Mirlin, Physics Reports 326, 259 (2000) 
^' M. Y. Azbel, Phys. Rev. B 27, 3852 (1983); M. Y. Azbel, Phys. 

Rev. B 28, 4106(1983). 

C. Basu, A. Mookerjee, A. K. Sen, and A. K. Thakur, J. Phys.: 
Condens. Matter 3, 9055 (1991). 
" J. B. Pendry, J. Phys. C: Solid State Phys. 20, 733 (1987). 

J. Bertolotti, S. Gottardo, and D. S. Wiersma, Phys. Rev. Lett 94, 
113903 (2005). 

W. Ren, F M. Xu, and, J. Wang, Nanotechnology 19, 435402 
(2008). 

C. W. J. Beenakker, Rev. Mod. Phys 69, 731 (1997). 



